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Abstract 

The two-body Coulomb scattering problem is solved using the standard complex scaling method. 
The explicit enforcement of the scattering boundary condition is avoided. Splitting of the scattering 
wave function based on the Coulomb modified plane wave is considered. This decomposition leads 
a three-dimensional Schrodinger equation with source term. Partial wave expansion is carried out 
and the asymptotic form of the solution is determined. This splitting does not lead to simplification 
of the scattering boundary condition if complex scaling is invoked. A new splitting carried out only 
on partial wave level is introduced and this method is proved to be very useful. The scattered part 
of the wave function tends to zero at large inter-particle distance. This property permits of easy 
numerical solution: the scattered part of the wave function can be expanded on bound-state type 
basis. The new method can be applied not only for pure Coulomb potential but in the presence of 
short range interaction too. 
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I. INTRODUCTION 



The method of complex scahng (CS) has been an excellent tool to calculate half life times 
of resonance states for a long time. The CS has been successfully applied in many areas of 
quantum physics |2| and it has been extended to collision processes very early on [3|, |4|. 
However, a drawback of the standard CS (or the uniform CS) has emerged immediately 
after the introduction of the method. For scattering problems the CS procedure can be 
applied only for short range potentials ^, This is indeed serious since the long range 
Coulomb interaction can not be neglected in majority of the problems of atomic and nuclear 
physics. Several modifications have been suggested but none of them has reached a 
widespread acceptance. After these initial applications the scattering aspects of the CS has 
been neglected. 

The turning point has been the work where it has been shown that scattering cal- 
culations with the exterior CS can be successfully performed for long range interactions. 
After this pioneering work the exterior CS method has been applied for variety of three- 
)ody Coulomb problems, even above the three-body breakup threshold, with great success 
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13| . The exterior CS method has proved to be one of the most successful numerical 



methods to deal with collision processes. However, recently the exterior CS method has 
been under scrutiny since in the method an artificial cutoff in some of the interaction is 
used. To solve this problem a modification of the original exterior CS method has been 
sugge^ted and checked ,n two-bod. calculations Q3 Extens.on to th.ee-body p.oblen, 



has been also sketched [H 

Recently it has been shown [17[ that the standard CS can be applied for scattering 
problems when a short range potential is added to the pure Coulomb interaction. The 
method is based on the two potential formalism. Similar approach has been suggested also 



m 



18|. In the present paper we rigorously develop a method which is equally good for pure 



Coulomb interaction and for the general case too (i.e. a short range potential is added to 
the Coulomb interaction). The new approach does not rely on the two potential formalism 
and dangerous cutoff will not be introduced. 

In the case of a two-body problem the wave function depends on the inter-particle coordi- 
nate r. The scattering solution of the Schrodinger equation with momentum k is denoted by 
ip'^(k,r). This wave function will be called three dimensional (3D) wave function. It is as- 
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sumed that the wave function satisfies appropriate scattering boundary conditions. The aim 
of the apphcation of any CS method is to introduce a new equation instead of the Schrodinger 
equation with simphfied boundary conditions. The expectation is that the solution of this 
new equation is square integrable therefore it can be approximated by bound-state type 
basis functions. In this way the exphcit use of the comphcated scattering asymptotic form 
of the wave function can be avoided and the numerical calculation can be simplified. 

In contrast to resonance state calculation in scattering problem the CS is not applied 
directly to the full wave function. First a splitting of the total wave function is carried out. 
The full scattering solution is searched in the form 

V^+(k,r) = 0o(k,r)+V^^'=+(k,r), (1) 

where 0o(k, r) is a known function. From the Schrodinger equation for the scattered part of 
the wave function the so called driven Schrodinger equation (or Schrodinger equation with 
source) 

{E-H)r'+ik,r) = Sik,r), (2) 

can be derived. The source term is given by S'(k, r) = {H — i?)0o(k, r). The Hamiltonian 
and energy are denoted by H and E respectively. 

We mention that the two-body Coulomb problem with source has been recently thor- 



oughly investigated in [19| . Complicated but exact solutions have been given for very general 
sources. Basis functions with proper two-body scatterin g as ymptotic have been generated 
from the exact solutions and used in the J-matrix method [20]. The driven Schrodinger equa- 



tion has been applied for realistic three-body scattering problems too [21|, |22| . However, in 
these works CS has not been applied and the complicated scattering boundary conditions 
have been implemented using either the finite element method or Sturmian expansion. 

The CS in scattering calculations means that the coordinate r in Eq. ([2]) is replaced by 
re'^ where < 6 < 7i/2. The boundary condition is simplified if after complex scaling the 
scattered part of the wave function goes to zero when the inter-particle distance tends to 
infinity. It is easy too see that this property is fulfilled if ip'^'^'^ {h, r) contains only outgoing 
spherical wave. In the paper this property will be investigated for different splittings of the 
full wave function. 

The organization of the paper is the following. In section |TT] the known expressions of the 
two-body Coulomb scattering are reviewed. The driven Schrodinger equation is introduced 



in section lllli The splitting of the total wave function in Eq. ([T]) is carried out on 3D level, 
however, the splitting can be carried out only on partial wave (p.w.) level. Different 3D 
and p.w. splittings will be considered and it will be shown that there are cases when the 
3D and p.w. splittings are not equivalent. The Coulomb modified plane wave (CMPW) 



plays a basic role in the recent surface integral formalism of the scattering theory [23|, |24 1 . 
The properties of the 3D splitting based on the CMPW will be investigated in section IIVI A 
useful p.w. splitting from the point of view of the CS will be introduced in section |Vl Finally 
numerical examples will be presented both for the pure Coulomb case and for a potential 
having short and long range parts. The conclusions will be given in section I VIII 

II. EXACT SOLUTIONS OF THE TWO-BODY COULOMB PROBLEM 

First we collect a few known expressions [25] for the two-body Coulomb scattering in order 
to fix the notations. As usual we take h = m = e = 1 {m is the reduced mass), the energy is 
E = k"^ /2 > and the Coulomb-potential reads 'jk/r, where 7 is the Sommerfeld-parameter. 
We consider the Schrodinger equation with pure two-body Coulomb interaction 

;Ar + ^-?V(k,r)=0, (3) 
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where Ar is the Laplace-operator. The Coulomb scattering state 

^'^(k, r) = e-'^^/2r(l + z7)e^'''"M(-z7, 1, ikr - zkr) (4) 



is a solution of ([3]). Here M{a,b,z) is the regular confluent hypergeometric function 26] . 
The partial wave expansion is given by the well known form 

00 

^+(k, r) = 5^(2/ + l)i:+{k, r)Pi{cos{}), (5) 
1=0 

where Pi{z) is the Legendre polynomial and is the angle between the vectors k and r. The 
full radial part il)^{k^r) is expressed with the help of the regular Coulomb function Fi{k,r) 

^+{k, r) = ^i' exp{tai)Fi{k, r). (6) 
kr 

The explicit formula reads 

= ^%^+2)^ e-^^''e-^'^{2ikryM{l + 1 - Z7, 2/ + 2, 2ikr) (7) 



and the Coulomb phase shift is defined by e^*'^' = r(/ + 1 + i'-f)/r{l + 1 — 27). The p.w. 
components tp^{k,r) satisfy the radial Schrodinger equation 



2r (ir^ 2r'^ 



^ r 2 



tpi{k,r) = 0. 



(8) 



The Coulomb scattering function can be split into so called incoming and scattered waves 



271]. Using the identity 7.2.2.9 in 28| we can write 



where 



and 



V'+l^r) =^i(k,r)+^,(k,r), 



iPi{k, r) = e^'^/^e^^^Ui-ij, 1, ikr - ikr), 



r(-«7) 



(9) 



(10) 



(11) 



The notation U{a,b,z) stands for the irregular confiuent hypergeometric function 26|]. In- 
terestingly not only '?/'^(k, r) but the functions ^ipiik, r) and ipsik, r) satisfy the 3D Scrodinger 
equation ([3]). 



The partial wave expansions of the incoming and scattered parts are given in 27| . Later 



we wi 



1 use them so we quote the main result of paper 



271]. We use a very similar notation 



as in 271] however we have rewritten the Whittaker function W in terms of U. 



The p.w. expansions of ijji{\<.,r) and ips{k,r) are given in the same form as (I5l) but the 
p.w. components now read 



i^iA^, r) = uji,i{k, r) + xi{k, r) 



(12) 



and 



V's,«(^5 r) = ujs,i{k, r) - xi{k, r). (13) 

We note that our definitions of ooi^i, ujg^i and xi ci-re constant times of the original ones 
The explicit expressions are the foUowings 

uji^i{k, r) = e-'^'e^^'^{-l)^+^{2ikr)^U{l + 1 - Z7, 2/ + 2, 2ikr), (14) 

u^ iik, r) = e'^'e^''''^^^'^{-iy+^{2ikr)^U{l + 1 + 27, 2/ + 2, -2ikr) (15) 

and 

(16) 



^^kr+,./2 r(2/ + 1) ^ (-1)"(Z7 - On ,,., 



2ikr {2ikry T{1 



n=0 



-2l)nn\ 



The equation 

ip^{k, r) = UJi,i{k, r) + ujs,i{k, r) (17) 

is also proved in We mention that the sphtting fll7p simply follows from ([7]) if the func- 
tion M{1 + 1 — i7, 2/ + 2, 2ikr) is rewritten in terms of the irregular confluent hypergeometric 



functions U using the equation 7.2.2.9 of 



28|. 



Although the functions '?/'i(k, r) and V^s(k, r) are solutions of the 3D Schrodinger equation 
([3]) surprisingly its partial wave components ■0jj(fc,r) and ips,i{k,r) do not satisfy the radial 
Schrodinger equation ^ (for details see [23|). This property will be proved to be very 
important for our new method. 

Using the asymptotic expansion 13.5.2 of 26|] we get the asymptotic form of the Coulomb- 



scattering wave function in the well known form 

1 



V^+(k,r) = e*''"(A;r - kr) 



kr 



^ikr—i-y ln(2fer) 



kr 



(18) 



where /c(cos(^9)) is the Coulomb scattering amplitude. The function e*'^''(A;r — kr)*"^ is called 
CMPW. 



III. DRIVEN SCRODINGER EQUATION 

The scattering solution of the Scrodinger equation is searched in the form ([1]). From 
the Scrodinger equation ([3]) with a simple rearrangement the following driven Schrodinger 
equation (or Scrodinger equation with source) 

C + ^^r-^V^^+(k,r) = 5(k,r), (19) 



2 2 r ^ 

can be derived for ^/'^^"''(k, r). The source term is given by 



5(k,r)=(-iAr + ^-^)0o(k,r). (20) 



2 r 2 

We mention that the driven Scrodinger equation (|T9l) has been studied in 19|. For quite 
general sources very complicated analytic solutions can be found [19!]. The aim of our paper 
is to derive an easy numerical method to solve (fT9ll and from the scattered part of the wave 
function deduce the scattering amplitude. 

The asymptotic form f|T8|) inspires the following choice for 0o(k, r) 

0o(k,r) = e*'^'"(A;r-kr)*^. (21) 



This splitting is based on the CMPW and it has been used in 23|, |2J| in order to derive the 
surface integral formalism of the scattering theory. Using Descartes-coordinates it is easy 
to derive a simple form for the source term 

7^^ 



5(k,r) 



e^'^^'ikr - kr) 
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(22) 



r{kr — kr) 

We may try to use the splitting based in the incoming Coulomb wave function i.e. we 
make the following choice 

0o(k,r)=^,(k,r) (23) 

instead of fl2T|) . In this case we get ^(k, r) = 0. This follows from the fact the function 
ipi{k,r) satisfies ([3]). In this case we do not get a driven Schrodinger equation, ip'^'^^(k,r) 
satisfies the original homogeneous equation ([3]) and ip^'^'^(k,r) = (k,r). 

If we want to derive the p.w. form of the driven Scrodinger equation (fT9l) we have to 
have the p.w. expansions of the source term 



^(k, r) = ^(2/ + l)Si{k, r)Pi{cos{^)) 



1=0 



and of 00 (k, r) 



6o(k,r) = 5^(2/ + l)(Po,i{k,r)Pi{cosm. 



Using the operator identity [29 1 



1=0 



A. 



(24) 



(25) 



(26) 



where is the square of the orbital angular momentum operator, we can derive the partial 
wave form of the driven Schrodinger equation f|T9l) 



1 _ + _ -fk 
2 2r dr'^ 2r'^ r 



ijr^{k,r) = Si{k,r), 



where 



iP''+{k, r) = ^(2/ + l)ijr~^{k, r)Pi{cos{^)). 



(27) 



(28) 



Later it will be proved to be very useful if we make the splitting of the scattering wave 
function not in the 3D form ([1]) but on the p.w. level. We take the p.w. component of the 
scattering wave function in the following form 



^+(/c,r) = 4)oAk,r)+^r+{k,r), 



(29) 



where (j)Q^i{k,r) is a fixed known function and il)l'^^{k,r) is considered as an unknown func- 
tion. From tlie partial wave Scfirodinger equation we get tlie following non-homogeneous 
differential equation 



1 dp 

1 r 

2 2r rfr2 



/(/ + !) 
2r2 



7A; 
r 



^r+(A;,r) = ^z(A;,r) 



where 



~Si{hr) 



^ r 2 



(30) 



(31) 



2r (ir^ 2r2 

If we take (f)o^i{k,r) identical to the partial wave component of 0o(k, r) i.e. (pQ^k^r) = 
(po,i{k,r) then the source terms Si{k,r) and Si{k,r) are identical if in equation ( !20!) the 
action of the Laplace-operator can be given in the form fl26l) . This replacement however is 
valid only for those functions which are finite at r = (see page 496 29|). 

In the case of splitting based on ( l23ll the function ipi{\<., r) is not finite at r = 0. In this 
circumstance the 3D splitting and the p.w. level splitting are different. We have already 
seen that the 3D splitting based on (!23|) does not lead to a driven Schrodinger equation. 
However if we make the following p.w. splitting 



^+(A;,r) = 7/'i,i(A;,r)+V^f^+(fc,r), 



(32) 



i.e we take (j)o^i{k,r) = ipi^k^r) then we get a driven radial Schrodinger equation. Direct 
calculation of (13T!) gives the following source term 

gifcr+77r/2 



Si{k,r) 



(33) 



2r2r(-i7)' 

Interestingly the source term is independent from /. The derivation of (!33|) is given in 
Appendix \M 



IV. PARTIAL WAVE EXPANSION AND ASYMPTOTIC FORMS 



The p.w. expansion of '?/'j(k, r) is given in 27|] and we have reviewed it earlier. We now 



determine the corresponding expansion of the CMPW. The p.w. expansion of the CMPW 
is written in the standard form 

00 

e^^^'ikr - kr)^^ = ^(2/ + l)r;(A;, r)Pi{cos{d)). (34) 
1=0 
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and the radial functions are given by the integral 

1 

ri^k,r) = -{kr-y^ e'^^'^l - xyWi{x)dx. (35) 
2 J -I 

A compact expression for the p.w. component of the CMPW can be given for arbitrary /. 
Using fl35l) and the integral 2.17.5.6 in |28|] we get 

n{k, r) = — t!2)i_(2A;r)'^e'*^'^ 2^^2(1 + 27, 1 + ^7; ^ + 2 + ^7, 1 + «7 - ^; -2tkr), (36) 

where (a)„ is the Pochhammer symbol. 

For the application of the complex scaling we have to know the asymptotic behavior 
of the scattered part of the wave function ipi'^^{k,r). Here we derive formulas valid at 
large r values. With the help of the expression 13.5.2 j26| we get the following asymptotic 
expansions valid at r — )■ 00 

^^^(k,r) ~ ^f^(-l)'+iy a^^^ (37) 

n=0 ^ > 

and 

' ' 2ikr ^ '•"(2^fcr)" ^ ^ 



The expansion coefficients are given by a'„ = (— 1)"(Z + 1 — «7)ri(— / — i^)n/n\ and a^n = 
(Z + 1 + 27)„(— / + i''))n/n\. These asymptotic forms show that with the help of the splitting 
( IT7|) the incoming and outgoing spherical waves are clearly separated in the p.w. Coulomb- 
scattering wave function. 

In order to derive an asymptotic expansion of Ti{k, r) we express (136!) in terms of Meijer's 



G function. Using 5.11.1(2) |30| we get 



Ti{k,r) = e*'(2fcr)*^(-l)'G2;3 I 2ikr 



(39) 



-Z7, — «7 
0, — 1 — ^7 — /, / — 

We are interested in the asymptotic behavior after the CS is carried out i.e. r is replaced 
by re*^ and < 9 < 71. We give the asymptotic expansion valid in this case. Considering 



the expression 6.5.32 



26| we can derive 



where the expansion coefficients satisfy the recursion 
4(n + l)dl^^ = 2\2n^ - n(2i7 -I- i-i\d\ - n{n -I - i^i - \){n^l - %i)d\_^ (41) 



0.02 



-0.02 




FIG. 1: The function ujifi{k,r) — To{k,r) and the next to leading order term of its asymptotic 
expansion ()44p are complex scaled using 9 = 0.1. Only the real parts are displayed. The momentum 
is /c = 3, the Sommerfeld-parameter 7 = 1/3 and / = 0. The solid line denotes the exact values 
and the dashed line corresponds to the asymptotic expansion (jH]). 

and (Iq = 1. 

Since the p.w. components are related to each other by the simple relation 

^l{k,r) = Ti{k,r) + ij^'^+{k,r) (42) 
and we have the splitting ( |T71) we can write 

iPr^ik^r) = [uJi,i{k,r)-Ti{k,r)]+UsAk,r). (43) 

We notice that the last term in fH3|) u)s^i{k,r) asymptotically contains only outgoing spher- 
ical wave (see Eq. fl38|) ) so the applicability of the CS is determined by the behavior of 
ujij{k^re^^)) — Ti{k^re^^)) at r — > oo. Fortunately the asymptotic expansions of the func- 
tions C0i^i{k,r) and Ti{k,r) are carried out using the same asymptotic sequence of functions 
{e'"^^^ {2krY^ [—1)^^^ / (2ikr)", n = 0, 1, 2 . . .} so we can simply add/subtract the asymptotic 
expansions as required [3l| . Using (1371) and (140|) we can write down the following asymptotic 
expansion 



[^i,z(^, 



re 



Ti{k,re 



ie 



2i/cre*^ 



{2kre'^Y<{~l)^+^Y^ 



n=0 



{2ikre'^Y 



(44) 



We realize that a - q — c^q = 0. This means that in leading or- 



Let's investigate (j3/^ v-^^Kr,^ fin^f _ 

der ijjf^~^[k,re^^) does not contain complex scaled incoming spherical wave. However in 
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higher orders ipi^'^ik, re*^) do contains complex scaled "generalized" incoming spherical wave 
(e~*'^''/r", n > 1). This means that the complex scaled scattered part of the wave function 
ip1'^^{k,re^^) does not tend zero as r — t- oo. 

This finding is demonstrated in Fig. [T] Both the left hand side and the right hand side 
of (1441) are displayed. From the asymptotic expansion only the next to leading order term 
is considered (the leading order term is zero). The real part of the function uJifi{k,re^^) — 
ro(fc,re*^) first starts to oscillate with decreasing order of amplitude however at larger r 
values the presence of the terms of the form ^-^^re''^ ^ ^^^le^n (dominate and the amplitude of 
the oscillation becomes larger and larger. 

We have got a very unfortunate result, if we use the 3D splitting based on the CMPW 
then the scattered part of the wave function asymptotically contains both incoming and 
outgoing spherical waves. This fact prevents the application of the complex scaling. 



V. COMPLEX SCALING AND SCATTERING STATES 



In the previous section we have established that the splitting of the wave function based 
on the CMPW i.e. the choice fl2T]) is useless from the point of view of CS. Now we turn 
to the splitting (!32l) which is carried out on the p.w. level. The scattered part of the wave 
function is given by 

2Pr+{k,r)=u;sAk,r)~Xi{k,r). (45) 

This equation follows from f|T2l) . f|T7j) and fl32|) . The asymptotic form fl38l) and the expres- 
sion (fT6|) for xi{k,r) shows that the scattered part of the wave function now contains only 
outgoing spherical wave and so the complex scaling can be safely applied. From f|T6|) and 
(!38l) we get in leading order 



{2kr) 



-17 



2ikr 



O 



1 



Jkr+j7r/2 (_l)'r(2/ + 1) 



2ikr 



\2ikr\J _ 
_r(2/ + l) 



o 



\2ikr\ 



(46) 



r(/ + 1 - Z7) 

Let's make a variable transformation and replace r with re*^ in the partial wave driven 
Scrodinger equation ( 1271) and furthermore introduce a new function with the definition 



i!f'g+{k,r) = e'^''/2^f+(fc, 



re 



(47) 
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where 9 is an arbitrary fixed real number. A simple calculation gives the following equation 

where the complex-scaled source term is defined by 

~Si,e{k,r) = e'^"^Si{k,re''). (49) 



The advantage of the complex-scaled driven Scrodinger equation fl48|) is that its solution 
behaves very simply asymptotically. If the scaling angle satisfies the condition < ^ < vr 
then from (HBl) it follows 

lim V5f^+(fc,r) = 0. (50) 

r— >oo ' 

From the asymptotic form fj46l) we can establish the following local representation of the 
partial wave Coulomb S-matrix 



^—ikre 



2tkre (fc, r)+ jv^ + ^ _ 



r — )• oo. (51) 



The local representation of the phase shift given in ij] is different from ( 15T|) since the 
splittings of the scattering wave function are distinct. 

The function 'ipl^^{k,r) is not regular at r = 0. However the validity of the limit 

limr'+i7/5f+(fc,r) = (52) 

r— i>0 

can be easily demonstrated. Details are given in Appendix [Bl In order to give simple bound- 
ary condition at r = we make the following transformation /ij*g^(A;, r) = r^^^ipi^g^lk, r). This 
transformation leads to regular function at r = 0. For the new function we get the following 
differential equation 

The price we pay for the simplification at r = is the appearance of first order derivative 
in the equation. 

From the earlier considerations presented it is obvious that the method based on the 
splitting (!32|) can be extended to case when a short range interaction Vs{r) is added to 
the pure Coulomb interaction. In this case the inhomogeneous differential equation flS5]) is 
replaced by 



- + e---— - e---- - e-- - K(re-) ) h^^ik, r) = r^^^Sj^,\k, r), (54) 

12 
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FIG. 2: The local representation of the phase shift for pure Coulomb potential (upper part) and 
for the general case when a short range potential is added to the Coulomb term (lower part). In 
the first case the exact solution (dashed black line) is also displayed. For the numerical solution the 
boundary condition is imposed at R=250 (red line) and R=1000 (green line). Detailed discussion 
is given in the text. 

where the new source term reads 

Sj^^{k,r) = Si,g{k,r) + e^3e/VM(fc,re*')K(re^'). (55) 

VI. NUMERICAL RESULTS 

The differential equations (153|1 or (154|) have to be solved with the boundary conditions 

ht:eHk,0) = (56) 

and 

lim h'i'g+{k,r) = 0. (57) 

r— >oo ' 
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In numerical calculations instead of ( 157|1 the boundary condition 

h'l^f+{k,R) = 0. (58) 

can be used. Here i? is a positive and otherwise arbitrary large number. The boundary 
condition f l58|) is of course an approximation and it has to be investigated how the result 
depends on R. The value of R should be in the asymptotic region where ( H6|) is satisfied. 

The finite element method is chosen as a numerical technique for the solution of Eqs . 
(15 3 p or The method and the basis functions used in any elements are described in 32 |. 



The same method was used also in [14!] . For the presented calculations equal 



y spaced finite 



32| is denoted by 



elements of length 1 is taken. The degree of the Lobatto shape functions 
N and the same N value is used at each elements. The 6 parameter of the CS was chosen 
to 0.1 radian. 

First the pure Coulomb case is considered i.e. the potential is given by 1/r. In this case 
the numerical result can be compared to the known analytical solution. The momentum 
was k = 3 and the considered orbital angular momentum was 1 = 0. The phase shift is 
calculated with the help of the local representation ( 15T|) . In this equation for ipi^g^{k,r) 
either the exact solution or the approximate one determined by the finite element method 
can be used. In the second example the potential 7.5r^exp(— r) is added to the previous 
pure Coulomb term. Exactly these two cases were studied in [14] where a different splitting 
of the wave function and the exterior complex scaling method was used. 

The upper part of Fig. |2] shows the results of the calculations carried out using pure 
Coulomb potential. In this case the exact solution (dashed black line) can be compared to 
the numerical ones. In the finite element method the boundary condition ( !58l) is imposed at 
two different R values {R = 250 and R = 1000). We note that for r > R the finite element 
solution is not defined. The boundary condition should be set at infinity (see (1571) ) but 
a finite R value is taken so it can be expected that the numerical solution is not accurate 
enough around the point where the boundary condition is set up. This can be clearly noticed 
in Fig. [21 If we choose R = 250 then there is an oscillation with large amplitude around 
r = 250. If the boundary condition is set up at a larger R value then the oscillatory region 
is pushed out around this value. In Fig. [2] the oscillatory region moved from r = 250 to 
r = 1000 simply changing the value of the parameter R from R = 250 to R = 1000. The 
effect of the boundary condition is noticeable. However, if this edge effect is not considered 
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r 



FIG. 3: The local representation of the phase shift for pure Coulomb potential. In the numerical 
solution the boundary condition is imposed at R=250. The exact solution and the numerical ones 
are displayed. In the numerical calculations the number of the Lobatto shape functions {N) is 
varied. 

then the local representation of the phase shift is practically constant on a huge region. This 
is a useful feature since it helps to determine a unique value of phase shift of the numerical 

n 

calculation. In contrast the local approximation of the phase shift in [1^] tends to the exact 
value by a persistent oscillation with decreasing order of amplitude. 

The lower part of Fig. |2] shows the results when the short range potential is added to 
the Coulomb interaction. This modification does not change the previous observations. The 
lower part of Fig. [2] clearly demonstrates each previous conclusions. The position of the 
boundary condition influences the value of the phase shift only around the point r = R. The 
phase shift calculated by the expression (I5T1) is practically independent from the value of r. 
We note that apart from the oscillatory region around r = 250 the two numerical solutions 
corresponding to the choices R = 250 and R = 1000 coincide in the region 50 < r < 250. In 
this region in Fig. [2] the red and green lines are indistinguishable on the used scale. 

The calculations displayed in Fig. [2] have been carried out using 50 Lobatto shape func- 
tions on each elements. We investigated the dependence of the local representation of the 
phase shift on the number of Lobatto functions used in the finite element method. The 
boundary condition (!58l) is set up at i? = 250. The results are depicted in Fig. [3l Apart 
from the region around r = 250 the exact phase shift is reproduced with three digits accu- 
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-0.1784 



0.1 0.2 0.3 0.4 

6 [rad] 

FIG. 4: The local representation of the phase shift for pure Coulomb potential as the function of the 
complex scaling parameter. Four different r values are considered and in the local representation 
of the phase shift the exact wave function is used. 

racy with = 50. The calculation with = 100 very well reproduces the exact solution 
(4 digits agreement is reached almost everywhere). 

It remains to check how the complex scaling parameter 6 influences the calculated phase 
shift. Figure H] displays the local representation of the phase shift as the function of the 
complex scaling parameter. Four different r values are used. In these calculation the exact 
wave function is used in fl5Tl) . If the value of r is in the asymptotic region (e.g. r = 550) then 
the calculated phase shift is independent from the value of the complex scaling parameter. 
For smaller r values it is advantageous to use larger 9 value to get better agreement with 
the exact phase shift. 

VII. SUMMARY 

We have rigorously shown that the two-body scattering problem of the pure Coulomb 
interaction can be solved using the standard complex scaling method. This is achieved 
without using any cutoff of the long range interaction. The intricate scattering boundary 
condition is greatly simplified and so the numerical solution can be plainly achieved. It 
is obvious that the suggested driven Scrodinger equation can be solved by the use of the 
exterior complex scaling method too. The advocated splitting of the total wave function 
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works for general circumstances. It can be applied not only for pure Coulomb force but short 
range interactions can be added to the Coulomb potential. It turned out that the splitting 
based on the Coulomb modified plane wave does not lead to simplification of the boundary 
condition from the point of view of the complex scaling. 



Appendix A 



In the case of the splitting fl32|) according to fl3T|) the source term reads 



Si{k,r) 



1 



2r dr' 



-J + 



/(/ + !) , 7fc k 
2r2 



(Al) 



Since Ui^i{k,r) is just the Coulomb Jost solution |25|, l27| the contribution from uji^i{k,r) to 
the source is zero. Introducing a new variable z = 2ikr and rewriting the summation in f|T6|) 
we have 



where 



Si{k,r) = 



P 2A;2 ^2 2k\l + 1) 2i-fk'^ 



-y+ . 



z — 



dz^ 

■yn/2 pZ/2 I 



r(i — i'f) z 



E 

n=0 



(/ + l)„(-0„(l)n^" 



Direct calculation gives 
Si{k,r) 



r(i-i7) z^ ^ (i-n) 

Rearranging the summation we get 



1 - i'y)n 

n{n + l) -/(/ + !) 



Siik,r) 



r(l-z7) z2 



n=0 



{l + lU-l)n 
(1 -^7)n 



+ 1) -/(/ + 1)] + 



+ (^7 - n - 1) 



+ l)ri+l(~On+l 



-71-1 



(A2) 
(A3) 

(A4) 



(A5) 



(1 - i-f)n+l 

Using the fact that the Pochhammer symbols satisfy the recursion (a)„+i = (a)„(a + n) we 
can show that the expression inside the square bracket is zero and so we have proved fl33l) . 



Appendix B 



According to ( 11311 in order to prove (1521) it is enough to show that 



lim r Us i{k, r) = lim r 

r— >0 ' r— >-0 



(Bl) 
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From the definition (fT6|) it follows 



r^o ' ^ (22A;)'+i r(/ + 1 - i7) ^ ^ 



Considering (|T5l) and the expression 7.2.2.3 28| the function r^^^Us,i{k,r) is a sum of three 
terms 

r'+^Us,i{k, r) = Ai{k, r) + Bi{k, r) + r), (B3) 

where 

A{k, r) = -^^(fe) '^^^r(2/"^lif(.7-^ log(-2zA:r)e-, (B4) 

and 

The following abbreviations are used 

= 12^^""'"'' 

and 

= + 1 + ^7 + - ^(n + 1) - ^(2/ + 2 + n). (B8) 

The digamma function is denoted by '^{z). From the expressions above it is clear that 
lim^-j^o ''^) = linir^o -B/(A;, r) = 0. The limit value of Ci{k,r) as r — )> is exactly 
the right hand side of (IB2p so we have proved fl52|) . 
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